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Abstract 

We study quasi-periodic tori under a normal-internal resonance, possibly with multiple eigen- 
values. Two non-degeneracy conditions play a role. The first of these generalizes invertibility 
of the Floquet matrix and prevents drift of the lower dimensional torus. The second con- 
dition involves a Kolmogorov-like variation of the internal frequencies and simultaneously 
versality of the Floquet matrix unfolding. We focus on the reversible setting, but our results 
carry over to the Hamiltonian and dissipative contexts. 

MSC-class: 37J40 

1 Introduction 

Persistence results for quasi-periodic motions were first proved for maximal tori in Hamilto- 
nian systems and became known as Kolmogorov-Arnord-Moser(KAM)theory. In [31] this was 
extended to lower dimensional tori and to other contexts like volume preserving and reversible 
systems. The role of the 'modifying terms' in terms of system parameters was clarified in |14ll24j 
and the Riissmann condition \13\ [35] allows to subsequently reduce the high number of param- 
eters to the bare minimum. 

These results yield what is called quasi-periodic (or normal linear) stability, i.e. families of 
invariant tori persist under sufficiently small perturbations when restricted to certain (measure- 
theoretically large) Cantor sets. The theorems in [13] make the crucial assumption that all 
eigenvalues of the matrix describing the normal linear behaviour be simple. This implies 
in particular that det ^1^0 (except for the disspative case and the high-dimensional volume 
preserving case, where this condition is explicitly added). Multiple resonances are admitted 
in [111 \T7\ |22] and the aim of the present paper is to admit zero eigenvalues without weakening 
the conclusion of quasi-periodic stability. 
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1.1 Setting and results 

We work on the phase space M = T" x E™ x M^p, where T" = (M/27rZ)" is the n-torus on 
which we use coordinates x = (xi, . . . , x„) (mod 27r), while on and M^^ we use respectively 
y = (yi, . . . , ym) and z = (zi, . . . , Z2p). In such coordinates a vector field on M takes the form 

x = f{x,y,z), y = g{x,y,z), z = h{x,y,z), 

or in vector field notation: 

X{x,y,z) = f{x,y,z)dx + g{x,y,z)dy + h{x,y,z)dz. (1.1) 

We assume that the vector field X depends analytically on all variables, including possible 
parameters which we suppress for the moment; referring to [14: \ 1241 [33] we note that our results 
remain valid when 'analyticity' is replaced by 'a sufficiently high degree of differentiability'. An 
invariant torus T of a vector field X is called parallel if a smooth conjugation exists of the 
restriction XIt with a constant vector field i; = w on T". The vector uj = {uji,uj2, . . . ,Xn) G 
is the (internal) frequency vector of T. The parallel torus is quasi-periodic when the frequencies 
are independent over the rationals. 

We are concerned with persistence of quasi-periodic tori under small perturbations, and to 
fix thoughts we concentrat^ on the reversible setting. To define reversibility we consider an 
involution (i.e. = I) 

G:M^M, {x,y,z)^{-x,y,Rz), (1.2) 
with R € GL(2p,M) a linear involution on M^^ such that 

dimFix(i2) = dim {z G M^p \ Rz = z] = p. 

The vector field X is then called G-reversible (or reversible for short) if 

G^{X) = -X. 

Using (jl.ip this reversibility condition takes the explicit form 

f{-x,y,Rz) = f{x,y,z), 
g{-x,y,Rz) = -g{x,y,z), 
h{—x,y,Rz) = —Rh{x,y,z), 

valid for all {x,y, z) € M. 

Following |12l [T3\ [T^ [2l] the vector field X is called integrable if it is equivariant with respect 
to the group action 

T«xM^M, (e,(x,y,z)) ^ (^ + x,y,z) 

of T" on M, or in other words, if the functions /, g and h in (II. ip are independent of the 
x-variable(s). Such an integrable vector field 

X{x, y, z) = f{y, z)d^ + g{y, z)dy + h{y, z)dz (1.3) 

^We give explicit formulations for reversible vector fields, but the results remain valid for e.g. dissipative, 
Hamiltonian or volume-preserving systems (vector fields and maps), where equivariance is also optional. 
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is reversible if 

f{y,Rz) = f{y,z), g{y, Rz) = -g{y, z) and h{y , Rz) = - Rh{y , z) (1.4) 

for all (y, z) G M™ x M^P; this implies g{y, z) = for all (y, z) eW^x Fix(i?). In case /i(0, 0) = 
thU n-torus To = X {0} x {0} is invariant under the flow of the vector field X. The normal 
linear part N{X) of (II. 4p at Tq is given by 

N{X){x,y,z) =ujd^ + nzd„ (1.5) 

with 

Lo = /(0, 0) and n = D^h{0, 0). 

We denote the subspace of infinitesimally reversible linear operators on M?p by gl_(2p;M) and 
by gl_|_(2|?;]R) the subspace of all i?-equivariant linear operators on R^^, i.e. 

Sl±(2p;M) ={n£ gl(2p;R) | QR = ±Rn}. 

In order to define the non-degeneracy of (II. 3p at the invariant torus Tq we consider the subspaces 

'^Hn = {^9, + ^^d, I e M", G Bl±(2p; M)} 

of the spaces of all G-reversible vector fields on M and of all G-equivariant vector 

fields, satisfying G^{X) = +X. For X G the adjoint operator 

adN{X):X — > X, Y ^ [N{X),Y] 

maps into X^"^; a similar statement is true for Xj^^ . 

Our interest concerns purely G-reversible vector fields, and G-reversible vector fields that are 
furthermore equivariant with respect to 

27r 2tti 

Fi : M — > M, [x,y,z)^{xi -,x^,y,zi,e i zjj). (1.6) 

Here zjj = Z2j-i + iz2j singles out two of the z- variables in a complexified form and zj = 
(zi, Z2, ■ ■ ■ , Z2j-2, Z2j+i, . . . , Z2p) Contains the remaining z-variables. To allow for a unified for- 
mulation of our results we define a reversing symmetry group S and a character (a group 
homomorphism) x ■ — *■ {il} as follows: 

(i) In the purely reversible case we set S := {Id, G} and xiG) ■= —1. 

(ii) In the equivariant-reversible case we define S as the group generated by G and Fi and 
define x by x{G) ■= —1 and xi^i) '■= 1- 

In both cases S is isomorphic to Z2 x Zi, the dihedral group of order 21. When I = 1 the 
generator Fi = Id of course is superfluous. For both cases we put 

X+ = {X G X \ E^{X) = X ior all E £ T,} 
X' = {X £ X \E^{X) = x[E)X iov all E gT.} 

^Often one has a whole family Ty — T" x {y} x {0} of invariant tori. While we are especially interested in 
bifurcations, the variable y will still act as a parameter, now unfolding the bifurcation scenario. 
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together with Xj^^ = n . Furthermore we let and B consist of the constant vector 
fields in X~^ and X~ ^ respectively and denote by 

O{no) = {Ad(A) • := AQQA~^\Ae GL+(2p;M)} 

the orbit under the adjoint action of GL-|_(2p;R) on gL(2p;M). 

Definition 1 (Broer, Huitema and Takens |14j ) The parametrise(^ vector field X\ with 
linearization N{X\){x,y,z) = uj{\)dx + ^{\)zdz is non-degenerate at \ = \q if 

BHT(i) keradAf(XAo) nS+ = {0}; 

BHT(ii) at \ = \q the mapping {u),VL) : W — > M" x 0l_(2p;R),A i-^ {uj{\),VL{\)) is transverse 
to MAo)} xO(0(Ao)). 

The two non-degeneracy conditions BHT(i) and BHT(ii) generalize the condition that a<lN{X\^^) 
has to be invertible, a requirement that lies at the basis of Mel'nikov's conditions ( p.7|) with 
|£| 7^ 0). One also speaks of BHT non-degeneracy. Compared to the formulation in |14j . § 8a2 
the requirement that 0(Ao) have only simple eigenvalues is dropped. The extension to multiple 
normal frequencies was developed in [I H 117 ^ l22] for invertible r2(Ao); we return to the original 
formulation of BHT(i). 

To formulate the strong non-resonance condition necessary for persistence of invariant tori we 
introduce for G gl_(2p;M) the normal frequency mapping a : Q\_{2p;R) — > M?P where the 
components of a{i}) are equal to the imaginary parts of the eigenvalues of $7 G 0l_ {2p; M). Higher 
multiplicities are taken into account by repeating each eigenvalue as many times as necessary. 

Definition 2 A pair (co, 0,) £ R*^ x gl_(2p; R) is said to satisfy a Diophantine condition if there 
exist constants t > n — 1 and 7 > such that 

\{k,u;) + {i,a{n))\>^\kr (1-7) 
for all ifc E \ {0} and i E Z^f with \£\<2. 

This condition is independent of the way in which we have ordered the components of a{0,); 
also, if {io, 0,) satisfies (II. 7p then the same is true for all {to, 0,) with ^1 E 0{0,). For each T G P 
we define the associated Diophantine subset 

F^ := {A E r I (lj(A), n{X)) satisfies ([LZl) } • 

When F is a small neighbourhood of some Aq E P where X is non-degenerate then F^ is nowhere 
dense but with large measure (provided that 7 is sufficiently small). 

Theorem 3 Let X E X~ be a family of T,-reversible integrable vector fields that is non- 
degenerate at Xq E p. Then there exists 70 > such that for aZZ < 7 < 70 the following 

^The role of the external parameter A occurring in Definition [1] can be (partially) taken by the internal 
parameter y. 
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is true. There exists a neighbourhood T of Xq, neighbourhoods y and Z of the origin in respec- 
tively andM?P, and a neighbourhood U of X in the compact-open topology on X~ such that 
for each Z ^lA one can find a mapping ^ ■.T'^ x y x Z x T — > M x P of the form 

^>(x, y, z, id, ^) = (^x + U{x, w, /i), y + V{x, y,uj,ii),z + W{x, y, z, w, /i), w + Aiiu, /i), + Miuj, /x) 

for which the following holds. 

(i) The mapping $ is T,-equivariant, real-analytic in the x-variable and normally affine in the 
y and z variables. 

(ii) The mapping <I> is C^-close to the identity and is a -diffeomorphism onto its image. 

(iii) The restriction of $ to the Cantor set T"' x {0} x {0} x of Diophantine X -invariant 
tori conjugates X to Z . The restiction of ^ to T" x y x Z x also preserves the normal 
linear behaviour to these invariant tori. 



In terms of [131114]. the conclusion of Theorem [3] expresses that the family X is quasi-periodically 
stable, i.e., structurally stable on a union of (Diophantine) quasi-periodic tori. This allows to con- 
dense Theorem [3 to the statement that non-degenerate E-reversible integrable vector fields are 
quasi-periodicall}0 stable. Quasi-periodic stability implies that for every small perturbation Z 
there exists a Z- invariant 'Cantor set' V C M x P which is a C°°-near- identity diffeomorphic 
image of the foliation T*^ x {0} x {0} x of n-tori. In the tori this diffeomorphism is an analytic 
conjugacy from X to Z, which also preserves the normal linear behaviour. 



1.2 Normal-internal resonances 

Resonances are at the core of the problems one has to solve when trying to prove quasi-periodic 
stability ~ persistence of elliptic invariant tori 

= T" X {y} X C Af := T" x M"" x R^p 

under small perturbation. The strong non-resonance conditions (jl.7p exclude in fact four types 
of resonances. An internal resonance 

{k, uj) =0 for some / /fc G Z" 

prevents the parallel flow on Ty to have a dense orbit whence the invariant torus is not a 
(minimal) dynamical object, but rather the union of closed invariant subtori. One cannot expect 
such an n-torus to persist, cf. \32\ I36j. for the same reason that a circle consisting of equilibria 
breaks up under perturbation (generically with only finitely many equilibria in the perturbed 
system). Such resonances are excluded by ()1.7p when taking £ = 0. 

For \£\ = 1 the inequalities ()1.7p constitute the first Mel'nikov condition, cf. [UEOlllI], and 
concern the normal-internal resonances 

{k,uj)=aj with fixed A; G Z" and i G {1, ... (1.8) 
*In [TT] one speaks of 'normal linear stability' instead. 
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Passing to co-rotating co-ordinates on yields this resonance with k = 0, ci. \10\ I16|. This is a 
2-step procedure. First k is brought into the form k = (ki , 0, . . . , 0) by means of a preliminary 
transformation 

N^N, {x,y,z) ^ {ax,y,z) (1.9) 

with a G SL{n, Z). For the second step we again write zj = {zi,Z2, ■ ■ ■ , Z2j-2, Z2j+i, • • • , Z2p), 
zji = {z2j-i, Z2j) and complexify zji = Z2j-i + iz2j. Then we perform a Van der Pol transfor- 
mation 

N-^N, {x,y,z)^ {x,y,zi,e"'''''zii). (1.10) 

The transformed vector field has a vanishing normal frequency aj = 0. Hence, already constant 
perturbations 

pd, = P2j-ldz,,_, + P2jdz,^ , (32j-l, P2j G K 

make the tori Ty move in a way that cannot be compensated on the linear level. Condition BHT(i) 
in Definition [T] prevents such perturbations whence Theorem [3] yields quasi-periodic stability, 
see also Corollary [6] in Section [3j An alternative to this condition is to take (generic) higher 
order terms of the unperturbed vector field into account. This typically results in bifurcation 
scenarios that turn out to be quasi-periodically stable (in an appropriate sense) as well, cf. O 

[IUl[2ni[2Il[MJ- 

The remaining possibility \i\ = 2 in (jl.7p excludes the normal-internal resonances 

{k, Lu) = a, lb aj with fixed A; G and i ^ j & {1, ■ ■ ■ ,p} (l-H) 

and 

{k,uj) = 2aj with fixed A; e Z" and j e {1, ... (1.12) 

For (jl.lip one can again achieve /c = in co-rotating co-ordinates, cf. 141], turning this normal- 
internal resonance into the normal resonance 

0/ai = ±aj, i / i G {1, . . . 

While now the invertibility of 0, does yield quasi-periodic stability of Ty , see [11^ [T7] and Corol- 
lary U] in Section [3l the normal behaviour still may be drastically affected. Using the y- variable 
as a parameter, e.g. the normal linear matrix 

/ 10 \ 

-10 

-1 
V 1 / 

in a conservative setting unfolds (or deforms) both to elliptic and to hyperbolic behaviour. Here 
it is the bifurcation scenario involving the surrounding tori of dimension n + 1 and n + 2 that 
can only be captured by taking higher order terms of the unperturbed vector field into account; 
quasi-periodic stability was achieved in [71 [9l [211 112] for the simplest conservative bifurcation 
scenarios. 

The remaining case p.l2p is meaningful only if not already implied by (jl.Sp , so assume that (jl.7p 
holds with \i\ < 1. Then we can still achieve A; = in co-rotating co-ordinates, but now on a 
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2-fold covering M — > defined as follows. The preliminary transformation (jl.Op brings the 
resonance vector k into the form k = 0, . . . , 0) with ki odd. The Van der Pol transformation 
is no longer a mapping from N to itself, but a covering mapping from M = T" x x R^^ 
onto defined by 

n : M — > N, {xi,x^,y,z) ^ (2xi, x*, y, z/, e'^'i^iz//). 

Here xi S and x^ = {x2, • • • , x^) G T""^. The deck group Z2 = {Id, F} of this 2-fold covering 
is generated by the involution 

F : M — > M, {xi,x^,y,zi,zii) ^ {xi - TT,x^,y,zi,-zii). (1.13) 

This means that 

no F = n. 

Note that this is the special case / = 2 of (jl.6p . the corresponding quasi-periodic stability is stated 
in Corollary [5] of Section [3l The resulting frequency halving (or quasi-periodic period-doubling) 
bifurcation scenarios are quasi-periodically stable in the dissipative [5] and Hamiltonian [21j 
settings and similarly reversible frequency-halving bifurcations may be expected to occur if 
appropriate non-degeneracy conditions on the higher order terms are fulfilled. 

In [irfl the second Mel'nikov condition ( (fTTl) with \i\ = 2) is avoided completely, i.e. also 
simultaneous normal-internal resonances (jl.lip and ()1.12p with differing k £ are allowed. 
The price to pay for this approach is that any control on the linear behaviour is completely lost. 
For instance, double eigenvalues itiai = ±102 generically unfold to a Krein collision, where an 
elliptic torus evolves a 4-dimensional normal direction of focus-focus type. Such changes cannot 
be captured without persistence of the (normal) linear behaviour. 

1.3 Contents and conclusions 

This paper fits in the framework of parametrised KAM theory [5l [71 [9l [T0l 113 ^ 121] that originates 
from Moser [31] ; in fact we present a generalization of [HI |l2l [H] , as well as of [TTl [22l [2l] . In 
the next section we explicitly work out two examples to which Theorem [3] applies. In Section [3] 
we elaborate the conditions of Theorem [3] and also formulate three corollaries; the Corollaries [5] 
and [6] make the novel results of this paper explicit. The proof of Theorem [3] is sketched in 
Section HI The necessary unfolding theory, which plays a key role, is deferred to the Appendix. 

Our approach allows for normal-internal resonances (jl.lip and ()1.12p (and possibly also (II. 8p ) 
with A; G Z" fixed. The ensuing deformations of the linear behaviour coming from the per- 
turbation are taken care of by considering a versal unfolding of the linear part i = flz of the 
unperturbed vector field, i.e., an unfolding that already contains all possible deformations. The 
necessary parameters are provided by y G M™'; the possibility that m > n distinguishes the 
reversible context from the Hamiltonian setting. An alternative is to let the system depend on 
external parameters A, where variation of (y, A) versally unfolds the linear part. 

The proof in Section H] is formulated in terms of filtered Lie algebras and therefore exceeds 
the reversible setting, carrying over to other contexts that can be formulated in these terms, 

^These papers consider Hamiltonian systems, but we expect the results to carry over to the reversible context. 
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notably the dissipative, volume preserving and Hamiltonian contexts; possibly combined with 
equivariance, cf. |1H I13j. In the Hamiltonian case this answers a conjecture formulated in [21] 
to the positive. For dissipative systems this has already been used in |5j when proving quasi- 
periodic stability of the frequency-halving bifurcation scenario. We expect that appropriate 
higher order terms in ()3.9p allow to obtain a similar result for reversible systems. 

The unfolding (jA.Sp recovers the result for the case p = 2 that was obtained in [25]. There 
a 4-dimensional reversible system with a codimension 2 singularity at the origin is studied by 
formal normal forms together with the persistence of the associated codimension 1 bifurca- 
tion phenomena. It would be interesting to investigate the persistence of the corresponding 
bifurcation scenario in the KAM setting. Note that an additional F-equivariance next to the 
G-reversibility would enforce the origin to be an equilibrium for the entire non-linear family, an 
assumption that is made in [25j . 

Acknowledgment. The second author received financial support by the European Com- 
munity's 6th Framework Programme, Marie Curie Intraeuropean Fellowship EC contract Ref. 
MEIF-CT-2005-515291, award Nr. MATH P00286. 



2 Applications 

We illustrate our results with two examples that explicitly show how our assumptions enter and 
what extra conclusions can be drawn. 

Example 1 (Quasi-periodic response solutions) To show how to check the appropriate 
assumptions we consider the simple example of a 1-parameter family of quasi-periodically forced 
oscillators 

z = ff,{t,z,z) = h^{t,ujt,z,z), (2.1) 

with a fixed frequency uj, for instance we take uj = ^(\/5 — 1) (the golden mean number). The 
forcing /i^ is 27r-periodic in the first two arguments. The search is for quasi-periodic response 
solutions with this same frequency vector {l,uj). 

Putting Zl = z, Z2 = z we can rewrite ()2.ip as an autonomous system 

Xl = 1 

±2 = U! 
Zl = Z2 

Z2 = h^{x,z) =h^{z) + h^{x,z) 

on X where we split hfj, into the average /i^ over x {z} and the oscillating part = 
h^ — h^. The integrable vector field given by 





= 1 


X2 


= uo 


Zl 


= Z2 


Z2 


= h^{z) 
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has invariant 2-tori for all G M with h^{zi, 0) = 0. These correspond to response solutions of 
the forced oscillator. 

Note that we allowed for to depend explicitly on Z2 whence Z2 i— > —Z2 is not a reversing 
symmetry. We impose the system to be reversible with respect to 

{xi,X2,Zi,Z2) ^ {-Xi,-X2,-Zi,Z2), 

in particular depends on zi only through zf and we concentrate on the invariant torus at 
z = 0. The dominant part 

N{x^)=d,,,+u;d,, + n{f,)zd, 

has the parameter-dependent 2x2 matrix 



1 



which is invertible whenever dih^{0) ^ 0. However, the non-degeneracy condition BHT(i) is also 
fulfilled if 9i/i^(0) = since the eigenvector to the resulting eigenvalue is not invariant under 
the involution 

■ -1 
1 



R 



From this we conclude that condition BHT(i) is always satisfied 
The non-degeneracy condition BHT(ii) is satisfied when 

d 

d/i 



5iV(0) / 0. (2.2) 



Thus, the system is BHT non-degenerate as soon as ()2.2p holds true. Therefore, given this by 
Theorem [3] (for an explicit formulation see Corollary ^ , if the oscillating part h is sufficiently 
small, the forced oscillator (12. ip has a response solution near z = 0, with linear behaviour 
changing where c?i/i^(0) passes through zero. □ 



Remarks. 

(i) Earlier for the existence of a response solution as an extra requirement the condition 
dihf,{0) / was needed dH [HI 131 [301 EI] . 

(ii) The stability change of the response solution as 3i/i^(0) passes through zero leads in the 
periodic case to additional periodic solutions bifurcating off from z = 0, cf. [551 [Ml 
We expect such bifurcations to carry over to the quasi-periodic case. 

We now return to the setting of the introduction where the normal-internal resonance ()1.12p led 
to a perturbation problem on a 2:1 covering space. The next example shows how the normally 
linear vector fields on the covering and the base-space relate to one another. 



9 



Example 2 (Multiple normal-internal resonance) On the phase space 



TV = X X = {x,y,z} 
we consider the normally linear vector field 

Y = 2dxi + ^dx^ + ^{^)zdz 

with 



/ 1 \ 

l + 1 

-^2 -1 -^1 

V -^2 i + m / 



where we think of the parameters = (w, //) G M'^ as been obtained from y G by localiza- 
tion (13. 7p . The eigenvalues ±1(1 + /ii) it \J—pi.2 of ri(/i) yield at /i = the normal frequency 
a = ibi that has two normal-internal resonances p.lip and p.l2p with the same k = (1, 0) G Z^. 
Complexifying both C/ = + 1^2 and = ^3 -|- i(^4 on the covering space 

N = M/(47rZ) X T X X = {6,6,^/,C} 
we have the covering mapping 

U:N^N, (6,6,^,0 ^ (eimod(27rZ),e2,r/,diag[e5'«i]C). 
This leads to the deck transformation 

F:N^N, (ei, ^2, r/, C) ^ (6-271, 6, r?,-C) (2.3) 
and the lifted vector field 

Y = (2)1% + (J2% + ^^(/^)Cf^C 

on N satisfying n*y = Y. In this setting ^ = ii, implying that cJi = 2 and the corresponding 
periods are Ti = 27r and Ti = vr, so Ti = 2Ti as should be expected. 

Regarding the Floquet matrices $7 and il. we have 

z = diag[^i6e5'«i]C + diag[e^^^i]C 
= diag[e5'«i] Qi^iC + j^c) 
= diag[e5'«i] (^ild+1^) C 
= diag[e^'«i] (^ild+Cl^ diagie'^'^^] z. 

Apparently 



n = diag[e5'«i] (^ild+fi) diag[e-^'^i] = ild+Q, 
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and the resulting family 



n{fi) = n{fi) - iid 



/ 
V 



of matrices is the LCU of 



-m 1 





-fJ'2 ^1 



/ 1 \ 

1 



V / 



\ 

1 

/ 



Every perturbation of y on can be lifted to a perturbation of y on iV that respects the 
deck transformation (|2.3|) and rescaling time we can always arrange xi = 2, i.e. that the first 
frequency equals 2. Applying Theorem [3] (for an explicit formulation see Corollary [5]) we may 
conclude that Y is quasi-periodically stable and this implies quasi-periodic stability of Y. □ 



It should be noted that such an application of kam Theory goes beyond the possibilities of 
im [T7t [22] . For n = 1 the full bifurcation scenario has been addressed in [6] and it would be 
interesting to develop the extension for periodic to quasi-periodic orbits. 



3 Main results 

In the perturbation problem we work on the phase space M = T" x x R^^' = {x, y, z] where 
we are dealing with a 'dominant part' 

X = uj, y = 0, z = fiz, or X = Lodx + ^zdz (3-1) 

in vector field notation. While it is always possible to translate a single given invariant torus to 
To = T" X {0} X {0}, it is an assumption on the system that this torus can be embedded in a 
whole family Tj, = T" x {y} x {0} of invariant tori parametrised by y. This can be equivalently 
stated as 

/i(y,0)=0 forahyeM"", (3.2) 

and the non-degeneracy condition BHT(i) in Definition [1] ensures that this assumption can be 
justified. For each e > the scaling operator 

p^:M^M, (x,y,z)^ (3.3) 

commutes with G and with the T"-action on M , and hence preserves reversibility and integra- 
bility. Using p.3p and the linearity of the push- forward (Dg)^^ of X under takes the 
form 

(P,),X(x,y,z) = V,{X{V;\x,y,z))) 

= f{^y, (^^z)dx + -g{ey, e^z)dy + -^h{ey, e^z)dz. 
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We can use (|1.4p to find that N{X) := lim (D^^j^X is the dominant part ()1.5p of X. The vector 

e— >0 

field N{X) = Lodx + ^zdz is again reversible and integrable; it is characterised by the frequency 
vector Lo = {loi, . . . ,ujn) G IR" which describes the fiow along the invariant tori Ty, and by the 
matrix 0, £ 0l(2p; M) which determines the linear fiow in the z-direction normal to the family of 
invariant tori. 

Since does not depend on the angular variable x E T" the vector field N{X) is in normal linear 
Floquet form. The Floquet matrix is infinitesimally reversible, satisfying 0,R = —RQ because 
of the reversibility of the vector field X. Observe that if ^ G C is an eigenvalue of 17 G gl_(2p; M) 
then so is —fi. Hence, the eigenvalues of Q can be grouped into complex quartets, conjugate 
purely imaginary pairs ±ia, symmetric real pairs and the eigenvalue zero with even algebraic 
multiplicity. 

3.1 Non-degeneracy conditions 

Since GL+(2p;M) is algebraic it follows that the orbit O(rio) is a smooth submanifold of 
gl_(2p;IR). The tangent space at rio to this orbit is given by 

rooO(Oo) = {ad{A) ■ Clo = AQo -QqAIAe gl+{2p;R)} = ad(Qo) (0l+(2p;M)) , (3.4) 

where we have used the fact that ad(A) • = — ad(r2) • A for all A,n G 0l(2p; M). An unfolding 
of 0,Q is a smooth mapping 

such that 0(0) = ^Iq. An unfolding is versal if it is transverse to 0{Qo) at /i = 0, which requires 
that s > codimO(Oo); a versal unfolding with the minimal number of parameters (i.e. with 
s equal to the codimension of O(rio) in 5l-(2p;IR)) is called miniversal. Using the Implicit 
Function Theorem it is easily seen that given a miniversal unfolding 0, : — > gl_ {2p; R) of 
^0 G Ql_{2p; R) we can write each J7 G Q\_{2p; R) near in the form Q = Ad{A) ■ 0(/u) for some 
{A, fi) G 0l_|_(2p;R) X R* close to (Id, 0) and depending smoothly on 17. In case all eigenvalues 
of are different from each other a miniversal unfolding amounts to simultaneously deforming 
the eigenvalues, see [H]. Our approach yields persistence results independent of the eigenvalue 
structure of (see [ID] for some other step towards such general persistence results). For more 
details on versal, miniversal (or universal) unfoldings we refer to [T| [2l [T9]. 

Property BHT(i) generalizes the invertibility condition required in the definition of non-degeneracy 
as it was formulated in |1H [12[ \T7\ [22] . What is really needed for the proofs is the invertibility 
of the linear operator 

a.dN{Xxo) -.6+ — >B- (3.5) 
and since dimFix(i?) = dimFix(— i?) this is fully captured by BHT(i). Computing 

adN{Xx,){m = -n{Xo)f3d, (3.6) 

shows that this certainly holds true if detr2(Ao) 7^ 0. However, the condition BHT(i) can still be 
satisfied if detr2(Ao) = 0, for example when ker(r2(Ao)) C Fix(— i?). The Floquet matrix r2(Ao) 
may have zero eigenvalues as long as the corresponding eigenvectors do not lie in B~^. 

Remarks. 
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(i) Up to now, the condition det l^o 7^ was one of the central assumptions for normal linear 
stability in the general dissipative context as well as in the volume preserving, symplectic 
and reversible contexts. Replacing this condition by BHT(i) allows to extend the known 
theorems to the singular case of eigenvalue zero. 

(ii) Property BHT(i) is persistent under small variation of A near Aq because of the upper-semi- 
continuity of the mapping A 1— > dimker J7(A). 

Property BHT(ii) means that locally the frequency vector w(A) varies diffeomorphically with A, 
while 'simultaneously' the local family A 1— > r2(A) is a versal unfolding of $^(Ao) in the sense 
of [U [2]. For earlier usage of this method in reversible kam Theory, see [T^ I17j . In the 
Appendix we develop an appropriate versal unfolding that depends linearly on A. 

When trying to answer the persistence problem for Ty it is convenient to focus on (a sufficiently 
small neighbourhood of) each of the invariant tori Ty [v G W^) separately, considering the 
label V £ M™ of the chosen torus as a parameter; formally this can be done by a localizing 
transformation, setting 

y = i^ + yioc and Xioc{x,yioc,z;u) := X{x,-u + yioc,z). (3.7) 

This way we get a parametrised family of reversible and integrable vector fields, still on the 
same state space M; in this localized situation we concentrate on the persistence in a small 
neighbourhood of the invariant torus Tq, corresponding to {yioc,z) = (0,0). For simplicity we 
absorb the additional parameter 1^ with the other parameters which we may have and we also 
drop the subscript 7oc'. 

The non-degeneracy condition BHT(ii) requires that s > n + codim O (J7(Ao)); in case all pa- 
rameters originate from a localization procedure this means that we should have m > n + 
codim O (r2(Ao)). Assume now that Xx is non-degenerate at Aq G M*, and let {luq,Q,q) := 
{uj{\o), il(Ao)). Using a re-parametrisation and a parameter-dependent linear transformation in 
the z-space we may assume that the parameter A takes the form A = {ui,fj,,fl) and belongs to a 
neighbourhood P of Aq := {ujq, 0, 0) in M" x x while the dominant part of the vector 

field reads 

N{X){x, y, z, w, fi, fj) = udx + Q{n)zdz, (3.8) 

where il. : M*^ — > gl_(2p;M) is a given miniversal unfolding of ^Iq. The /i-part of the parameter 
does not appear in this expression for the (unperturbed) vector field X; although it might appear 
explicitly in the perturbations it turns out that /i plays no role at all in the further analysis. 
Therefore we suppress jl and just keep the essential parameters (w, n) and set s = n -|- c, with 
c = codim 0(17(0)). The question of a particular choice for the miniversal unfolding ^{fi) 
appearing in ()3.8p is addressed in the Appendix. 

3.2 Consequences 

We are given a family of integrable vector fields 

X{x, y, z, io, ^) = [uj + f{y, z, lo, fi)] + g{y, z, lo, fi) dy + [Q{fi)z + h{y, z, cj, /i)] 9^ (3.9) 
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on the product M x P of phase space M = x M™ x M?^ and parameter space P C M'^ = M" x 
with reversing symmetry group T, generated by (jl.2|) and (|1.6p . For I = 1 the latter is just the 
identity, but for / > 2 the composition 

27r 

Hi:=FioG:M-^M, {x,y, z) ^ {— - xi, x,,y, SiRz), (3.10) 

is another reversing symmetry and one may also characterise the vector fields in X~ as being 
reversible with respect to the two mappings G and Hi. Note that in this characterisation Hi 
may be replaced by o G for any i relative prime to /. 

The coefficient functions /, g and h entering X are higher order terms in z, satisfying f{y,0,uj,fj,) = 
g{y, 0, uj, fi) = h{y, 0, /i) = Dzh{y, 0, uj,^) =Q for all y G R™, (u, ^) G P. Within X~ we con- 
sider perturbations Z of X and write 



Z{x,y, z,u!, n) = UJ + f{x,y, z,iLi, fi) dx + g{x,y, z,tLi, fi) dy + Q{fi)z + h{x,y, z,uj, /j,) 



d 



here the coefficient functions /, g and h may contain lower order terms but are close to /, g and h, 
respectively. In this situation Theorem [3] allows to conjugate Z to X as far as Diophantine tori 
are concerned. 

The condition that $ be a full conjugation from X to Z means that <5*(X) = Z, or equivalently 
(<I>~^)^ (Z) = X. What is actually proved is the existence of a local diffeomorphism ^ such that 

{^-^)jZ){x,y,x,i^,fi)=N{X){x,y,z,uj,i^)+0{\y\,\z\)dx + 0{\yl\z\^^)^^ 

(3.11) 

for all {Ld,n) G Py which are sufficiently close to {ljo,0). The property ()3.1ip implies that for 
all parameter values {to, n) in the indicated Cantor set the X-invariant torus T" x {0} x {0} is 
mapped by $ into a Z-invariant torus on which the Z-Qow is conjugate to the constant flow cudx 
on This means that a Cantor subset of large measure of the family T" x {0} x {0} x P of 
X-invariant tori survives the perturbation to Z. The preservation of the normal linear behaviour 
means that the normal linear vector flelds X(X) and N[Z) along two corresponding invariant 
tori are conjugated by the derivative of the C°°-near-identity diffeomorphism. 

In comparison to earlier results on persistence of lower-dimensional tori the condition that all 
eigenvalues be simple is dropped in Theorem [3] and the condition det $7(0) 7^ is weakened 
to BHT(i). Indeed, we have the following corollary. 

Corollary 4 (Ciocci [17] . Broer, Hoo and Naudot [11]) Let the family X G X~ of G- 
reversible integrable vector fields satisfy the nan- degeneracy condition BHT(ii) at Xq = {ujo,0) G 
P, with 0(0) invertible. Then X is quasi-periodically stable. 

Next to the above purely reversible case I = 1 also the case Z = 2 of a reversing symmetry group 
E = {Id, F, G, H} merits an explicit formulation. Here H = H2 is given by (j3.1Up and yields 

f{y,SRz) = f{y,z), g{y, SRz) = -g{y, z) and h{y, SRz) = -SRh{y, z) (3.12) 

for integrable vector fields where S{zj,zij) = {zi,—zjj). From (??) it follows that f,g are 
even in zu, while h is odd in zij. Moreover, (jl.4p and (|3.12p imply that g{y,z) = for all 
{y, z) G M"" X Fix(P) and also for all (y, z) G M"" x Fix(5P). 
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Corollary 5 Let X G be a family of G-reversible F-equivariant integrable vector fields 
that satisfies the non- degeneracy condition BHT(ii) at Aq = {^^o,0) € P, with Q(0) invertible 
on Fix(5). Then X is quasi-periodically stable. 

Again we allow for multiple eigenvalues, in particular the eigenvalue may have multiplicity 
larger than two. A similar statement holds in case of equivariance with respect to (II. 6p instead 

of (HH. 

In the covering setting of Section 1, we observe that the lift of an integrable vector field again is 
integrable. In fact, if X is the lift to M of an integrable vector field X on N, then 11* (X) = X 
and F^X = X. In case the second Mel'nikov condition is violated by a resonance (|1.12|) we can 
apply Corollary [5] on a 2:1 covering space. In Example [2] of Section [2] we do this for a double 
normal-internal resonance with fixed resonance vector k ^1?. 

Corollary 6 Let X E X~ be a family of G-reversible integrable vector fields that satisfies the 
non- degeneracy condition BHT(ii) at Aq = (wojO) S P. //kerr2(0) is contained in Fix(— i?) then 
X is quasi-periodically stable. 

If kerr2(0) C Fix(-|-i?) we generically expect a quasi-periodic centre-saddle bifurcation to take 
place, cf. [20] . Here violations of the first Mel'nikov condition prevents persistence of the cor- 
responding tori if not embedded in an appropriate bifurcation scenario. The scaling (|3.3|) also 
can be applied to non-integrable systems, making the non-integrable higher order terms a small 
perturbation. It is then not automatic that the resulting dominant part is in Floquet form, this 
is a necessary extra requirement that can be thought of as generalization of integrability under 
which quasi-periodic stability can still be achieved. For a more thorough discussion of these 
questions see 

4 Sketch of proof 

The proof of Theorem [3] follows [11[TT1[13] almost verbatim (see also [12l [T71 [22] ) and here we just 
concentrate on the novel aspects. The quite universal set-up of [lll[3l] is based on a Lie algebra 
approach, using a standard Newtonian linearization procedure. The conjugation $ between the 
integrable and the perturbed family is produced as the limit of an infinite iteration process. The 
central ingredient of the proof is the solution of the linearized problem, the so-called homological 
equation. The structure at hand, that is, the reversible symmetry group S, is phrased in terms 
of the Lie algebras Xj^^ and and therefore automatically preserved. Here we content 
ourselves showing how the non-degeneracy conditions BHT(i) and BHT(ii) enter when solving the 
homological equation. 

At each iteration step we look for a transformation (^, r/, a, i^) {x,y,z,uj,fi) with lo = 
a + Ai((T, I/) and n = v + A2{a, v) independent from the variables (^, r/, C,) so that the projection 
to the parameter space P is preserved. The transformation in the variables is generated by a 
S-equivariant vector field ^' G X^ that we write as 

= ud^ + vdy + wd,. 

The homological equation reads 

SidN{X){^) = L + N (4.1) 
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with 

L^A^,r],C) = {Z -X}iin,d and A^^,^(e, 77, C) = Ai(ct, z.)^^ + 1^(A2(ct, i/))C5c 
and determines the unknown U, V, W, Ai and A2 according to 

U^a = Ai + 7 

V^a + V^n{u)C = g + Qr^ri + gcC (4.2) 
W^a+[n{v)C,W] = n{A2)C + h + h^rj + h^C- 

Here U,V,W, f ,g, h and their derivatives depend on 0, 0, a, i^). Moreover, greek subscripts 
denote derivatives, while U^a = Ti^^iU^^Oj and similarly for V and W . These linear equations 
are solved by suitably truncated Fourier series. Note that the left hand side of ()4.2p consists of 
the components of the vector field ad A^(Xo-,,y)(^), where 

For a given Z (and hence L), the goal is to find ^ G '^und ^ ^ keradA^(X)^ C X^^^ so 
that the homological equation ()4.ip is satisfied. Here ^ = X^^ n X^ denotes the intersection 
set of the Taylor and Fourier truncations of vector fields in X^ . 

We make the ansatz 

V{i,7],(:,a,v) = Vo + Vii] + V2(: and iy(e, r/, C, a, z^) = Wq + PFit? + W2C (4.3) 

for the unknown ^, where Vj and Wj (j = 0, 1, 2) depend on ^ and on the multiparameter (cr, v). 
Fourier expanding in ^ and comparing coefficients in (|4.2p yields the following equations for an 
explicit (formal) construction of To avoid clumsy notation we suppress the dependence on 
{a,v). 

For A; / 0, equation (14. 2p implies 

i{k,a)Uk = Ik (4.4) 
i(fc,cr)^o,fc = (4.5) 
\{k,a)Vi,k = {g^)k (4.6) 
y2,fc[i(A:,^)Id+J^(z^)] = (5c)fc (4.7) 
[i(A;,cT)Id-0(i/)]W^o,fc = hu (4.8) 
[i(A:,a)Id-17(i/)]TFi,fc = (/ijfe (4.9) 
[i(A;,a)Id-adf^(i/)]VF2,fc = (hi^)k (4.10) 



and, similarly, for fc = 



-Ai = /o (4.11) 

V2,oJ^(i^) = (9c)o (4.12) 

-fi(z.)T^o,o = /^o (4.13) 

-^{v)Wi^^ = (hr,)o (4.14) 

adn{iy)W2,k - ^(^2) = {h()o. (4.15) 
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On the one hand, it is clear by the Diophantine conditions that for k ^ none of the coefficients 
at the right hand sides of (j4.4|) - (|4.9p is in the kernel, i.e. none of the eigenvalues i{k,a), i{k, a)zizXj, 
with Xj eigenvalue of $7(z^) are zero. For k = 0, the equations (I4.12p - (l4.14p are solvable by the 
non-degeneracy condition BHT(i) since the right hand sides lie in B~ . The so-called solvability 
condition ()4.11|) determines the 3^-component 

of N in (14. ip . Turning our attention to equation (I4.10p . we see that it admits the solution 

W2,k = ma)ld-adn{iy)]-\h^)k 
if and only if the operator [i(A;, a) Id — ad ^{i^)] is invertible, which boils down to the condition 

i{k, a) / Xj - Xi 

on the spectrum of adr2(z^), where Xj is an eigenvalue of 17(1/). This inequality is the second 
Mel'nikov condition and again guaranteed by the Diophantine conditions. For k = the splitting 

im(ad+(rJo)) eker(ad_(J7^)) = gL(2p,M), (4.16) 

detailed in the Appendix lies at the basis of solving equation (j4.15p . Indeed, the non-degeneracy 
condition BHT(ii) guarantees that we may choose the LCU for fi. Using the Implicit Function 
Theorem and the fact that O (by construction) is an isomorphism between parameter spaces, it 
follows that ()4.15p admits the solution 

A2{a,u)=n-^(^-7T(h^,o + adn{i^)W2,o)), (4.17) 

where the mapping vr denotes the projection of gl_(n,M) onto the subspace ker(ad_(r2^)) ac- 
cording to the splitting ()4.16p . Compare with [I7j . Lemma 8.1. 

A Unfolding reversible linear matrices 

Let fio € gl_(2p;M) be given; the aim of this appendix is to summarize some results from 
[IZlEHlETj which allow to describe a miniversal unfolding of JIq) and to work out the details for 
two particular cases. 

Let r^o = Sq+Mq be the Jordan-Chevalley decomposition of f^o into commuting semisimple and 
nilpotent parts. The uniqueness of this decomposition implies that both Sq and A/q belong to 
0l_(2p;M). Also 

ker ad{Qo) = ker ad(5o) n ker ad(7Vo), (A.l) 
as easily follows from the fact that Sq and A/q commute, and as shown in |17^ I27j furthermore 

ker ad(O^) = ker ad(cSo) n ker ad{M^). (A.2) 

We know from (j3.4p that T^gO{Qo) = im (ad+(ilo))) while a classical result from linear algebra 
shows that the subspace ker (ad_(f]Q )) of gl_(2p;M) forms a complement of the tangent space 
Tn^p{no) to the orbit through ^q. Finally, ker {ad-{nl)) = ker (ad_(cSo)) n ker (ad_(A/'o^)) by 
()A.2p . and hence we obtain the following result. 
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Theorem 7 Let ^Iq G gl_(2p;M) be given, and let CIq = Sq + A/q be the Jordan- Chev alley 
decomposition of Qq . Then 

O : ker (ad_(cSo)) nker (ad_(A/cf)) — > Ql_{2p;R), A ^ + A, (A.3) 

forms a miniversal unfolding of £ gL {2p; R) . 

The unfolding is in the centrahzer of Sq. In the present context of linear systems one calls 
such an unfolding a linear centrahzer unfolding (lcu for short). Also note that r2(/u) — l^o is 
linear in the unfolding parameters. For the convenience of the reader we now explicitly work 
out a linear centrahzer unfolding (|A.3p for three particular choices of {Qo,R). 



A.l Unfolding multiple non-zero normal frequencies 



For our first example we assume that l^o G gL(2p;IR) has a 1 : 1 : • • • : 1 resonance (or p-fold 
resonance), meaning that 0,q has a pair of purely imaginary eigenvalues, say iti, with algebraic 
multiplicity p; we furthermore assume that we are in the generic situation, with geometric 
multiplicity 1. The subspaces ker(A/Q) (1 < j < p) form a strictly increasing sequence of 
subspaces invariant under Sq and R, with dimker(A/o) — dimker(A/o~^) = 2. With respect to a 
conveniently chosen basis {u^ } of the linear matrices Qq and R have 

the matrix form 



( J2 

O2 



V 02 



J2 
J2 



02 
J2 



02 \ 

02 

J2 
J2 / 



R 



( R2 O2 O2 
O2 R2 O2 



V O2 ... 



with 



1 
-1 



O2 








and Ro 



... O2 \ 

■•• O2 

■•• O2 

O2 R2 J 



From this one can compute an LCU of 0,q as follows. 



Lemma 8 Fix an A G ker (ad_(5o))nker (ad_(A/Q")). Then there exist constants /ii, /i2, 
M such that if we set 



(A.4) 



(A.5) 



A, :=A-Y, (AAo^ 



(1<J<P), 



i=l 



then Aj(Up^i) = {0} for 1 < j < p and < i < j — 1. In the particular case that j = p we have 

p 



A 



i=l 



JlN. «— 1 



(A.6) 
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Combining ()A.4p and ()A.6P an LCU of ^Iq takes the explicit form 



/ /ilJ2 O2 O2 
^2J2 ^lJ2 O2 



\ AipJ2 



... O2 \ 

■•• O2 

■•• O2 

fl2^2 ^lJ2 / 



(A.7) 



with unfolding parameters /ii, . . . G M. This construction invariably leads to the same LCU, 
we therefore speak from now on of the LCU. In case all eigenvalues of ilo G Ql_{2p;M) are 
purely imaginary, non-zero and with geometric multiplicity 1, the LCU of can be obtained by 
considering the different pairs of eigenvalues ±iaj, multiplying ()A.7p with aj (using for each j the 
appropriate dimension and a new set of parameters), and juxtaposing the obtained unfoldings 
as blocks along the diagonal. 



A. 2 Unfolding multiple eigenvalue zero 

For our second and third example we assume that 0,q € 0l_(2p;IR) has as an eigenvalue with 
geometric multiplicity 1 and algebraic multiplicity 2p; then Sq = 0, A/q = A/q 7^ for 
1 < J < 2p, and Mq^ = 0. The subspaces keT{J\f^), 1 < j < 2p are invariant under R; they 
form a strictly increasing sequence, with dimker(A/'Q) — dimker(A/'g"^) = 1. With respect to 
a conveniently chosen basis = A/q is a classical nilpotent Jordan matrix with I's above the 
diagonal. The matrix form of R depends on whether ker(A/'o) C Fix(i?), in which case R has 
the same matrix form as in (|A.4p . or ker(A/'o) C Fix{—R), whence the matrix form of R equals 
minus the expression in ()A.4|1 . 

To determine the LCU of JIq we first consider some A G ker(ad(7V^)); one easily shows that A 



can be written as A = X]j=i {■^0)'' 
the further condition that A E qI_ {2p; ] 
we obtain then the following LCU 



with some constants I'j G M (1 < j < 2p). Imposing 
gives Vj = for j odd; setting fj.j := i>2j for 1 < j < p 



Hence ilg has co-dimension p and the LCU is given by 



/ 1 
10 
10 



\ 






1 



+ 



= 






...,fip) G MP. 


by 















\ 


^1 












Ml 







^2 





Ml 






















fl2 Ml / 



(A. 
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alternating diagonals with unfolding parameters fij and diagonals with 0. Note that we may 
alternatively fix R to be of the form (|A.4p and obtain the two cases by taking (|A.8p and its 
transpose, with 0,q = A/q having its I's below the diagonal. 

In case the condition dimker(r2o) = 1 on the geometric multiplicity of the zero eigenvalue 
is dropped the unfolding changes drastically and requires more parameters, i.e., has higher 
codimension. The same is true for our first example (non-zero normal frequencies). Further 
information on these cases can be found in |23| . 
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